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Abstract 

a 

o 

Q | We discuss the low-temperature behavior of the electronic self-energy in 

the vicinity of a ferromagnetic instability in two dimensions within the two- 
C*~> ' particle self-consistent approximation, functional renormalization group and 

r» ' 

t> 
o 

t^J- ' at T > 0, the self-energy has a non- Fermi liquid form at low energies |w| < Ag 

O ' 

near the Fermi level, where Ao is the ground-state spin splitting. The spec- 

c3 ; 

tral function at temperatures T < Ao has a two-peak structure with finite 
spectral weight at the Fermi level. The simultaneous inclusion of self-energy 

o' 

Q ' and vertex corrections shows that the above results remain qualitatively un- 

changed down to very low temperatures T -C A . It is argued, that this form 
of the spectral functions implies the quasi-splitting of the Fermi surface in the 
paramagnetic phase in the presence of strong ferromagnetic fluctuations. 
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I. INTRODUCTION 



Non-Fermi-liquid behavior of correlated low- dimensional electron systems has attracted 
much attention during the last decade. This behavior is usually connected with the viola- 
tion of the quasiparticle (qp) concept in some energy window around the Fermi level. A 
prominent example is the pseudogap phenomenon observed in underdoped high-T c com- 
pounds [1]. Cuprate superconductors are however not the only materials example which 
show a strong suppression of the low energy spectral weight due to correlation effects. A 
qualitatively similar behavior on parts of the Fermi surfaces was also observed recently in 
the unconventional superconductor Sr 2 Ru04 in an intermediate temperature range [2] and 
the layered manganite compound Lai +x Sr2- x Mn 2 07 [3]. 

One possible viewpoint for the origin of the pseudogap in high-T c compounds is to 
relate it to precursors of antiferromagnetism [4-6]. The fluctuation exchange (FLEX) [7,8], 
two-particle self-consistent (TPSC) [9], dynamical cluster approximation [10], and most 
recently the functional renormalization-group (fRG) technique [11,12] have demonstrated a 
strong anisotropy of spectral properties around the Fermi surface (FS) in the 2D Hubbard 
model with a possible violation of the qp concept on parts of the FS. The large incoherent 
contributions to the electronic spectrum of low- dimensional metallic antiferromagnets (e.g. 
an anomalously large scattering rate) were also discussed in Ref. [13]. These phenomena 
result in the formation of a two- [7-9,12] or three-peak [11] structure of the spectral function 
in the vicinity of the antiferromagnetic (AFM) instability. With decreasing temperature, 
these pseudogap features are expected to evolve continuously towards a ground state spectral 
function with an AFM energy gap at the Fermi level. 

While many results have been obtained for the electronic properties in the vicinity of 
an AFM state and some results exist for itinerant ferromagnets [14,15], surprisingly much 
less is known about the evolution of the quasiparticle properties in the paramagnetic phase 
near the ferromagnetic (FM) instability. The paramagnon theory [16-19] focused mostly on 
the description of the magnetic properties, ignoring the renormalization of the one-particle 
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Green function. Nevertheless, already Doniach and Engelsberg [20] showed that in three 
dimensions the qp weight vanishes logarithmically on approaching a FM zero-temperature 
quantum phase transition (QPT). For two-dimensional (2D) systems, an e 2 ^ 3 energy depen- 
dence of the self-energy at the quantum critical point (QCP) can be inferred from similar 
calculations in the context of gauge field theories [21], the phase separation problem [22], 
and the Pomeranchuk instability [23]. The latter two instabilities arise also in the zero- 
momentum transfer particle-hole channel and therefore are expected to have properties 
which are similar to those in the vicinity of the ferromagnetic QPT, at least in lowest order 
perturbation theory with respect to the fermion-boson (charge or spin modes) coupling. The 
£ 2/3 f re q UenC y dependence of the self-energy implies the vanishing of the qp weight at the 
Fermi level and therefore invalidates the qp concept. This may raise doubts about the va- 
lidity of the applied scheme, since the above-mentioned calculations did not consider both, 
self-energy and vertex corrections. However, the calculations by Altshuler et al. [24] within 
the gauge field theory context showed that the e 2//3 dependence of the self-energy remains 
valid also in higher orders of perturbation theory. 

The breakdown of the qp concept at the QCP may be even more apparent at finite 
temperatures. For fermions interacting with a gauge field it was shown that for the case of 
gapless (although diffusive) bosonic excitations the imaginary part of the self-energy in a 
non-self-consistent calculation is divergent at the Fermi level at T > as a consequence of 
the diverging static spin susceptibility x(0, 0) [21]. This divergence should necessarily have 
certain consequences for the zero-momentum particle-hole instabilities of fermion systems 
with short-range interactions. Although the magnetic correlation length £ of 2D systems 
with short-range interaction is finite at finite T and x(0, 0) oc £ 2 , three temperature regimes 
should be distinguished [25] (see Fig. 1): (i) the quantum disordered regime with a disor- 
dered ground state and almost temperature-independent correlation length, (ii) the quantum 
critical (QC) regime with £ ~ 1/T, and (iii) the renormalized-classical (RC) regime above 
an ordered ground state with an exponentially large correlation length at low T. The di- 
vergence of the imaginary part of the self-energy for £ — > oo may lead to especially strong 
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effects in the RC regime, where the inverse correlation length is almost negligibly small at 
low T. Indeed, for the AFM instability this divergence results in the formation of a pseu- 
dogap structure in the spectral function [9]. Although similar properties in the vicinity of 
a FM instability were discussed quite recently [26,27], the behavior in this case is far less 
from clear, since the suppression of the spectral weight at the Fermi level itself weakens 
the tendency to ferromagnetic order, which therefore implies the necessity to account for 
self-energy and vertex corrections selfconsistently. 

The recently found non-analytic contributions to the spin susceptibility in second-order 
perturbation theory with respect to the electron-electron interaction [28] have led to infer 
the possible absence of a second-order FM QPT, which is either replaced by a first-order 
transition to ferromagnetism or a second-order QPT into an incommensurate phase with a 
finite ordering wavevector Q, which then continuously decreases towards Q = on moving 
away from the QPT into the ordered state. However, these corrections are not expected, at 
least in the weak-coupling regime, to remove the renormalized classical temperature regime 
entirely Indeed, the corresponding characteristic temperature scale T x ~(£//47rt) 2 t, below 
which these corrections become important [28], is quadratic in the interaction and therefore 
small in the weak-coupling regime in comparison with both the bandwidth and the crossover 
temperature T* ~ U into the renormalized classical regime not too close to the QPT. 

It therefore appears demanding to investigate the finite-temperature behavior of the 
self-energy in the vicinity of the FM instability and to compare non-self-consistent and self- 
consistent techniques. For non-self-consistent calculations, the TPSC as well as the recently 
proposed fRG approaches on a patched FS [29-32] can be used. The latter approaches 
have the advantage that they do not select certain types of electronic scattering processes, 
but consider them all on equal footing. The self-consistent treatment of self-energy and 
vertex corrections is per se a rather difficult task. Most consistently this can be done in 
a parquet-type analysis [34] (see also Refs. [35,36]), which is however necessarily rather 
involved for numerical studies in dimensions d > 1. In the fRG technique the back influence 
of the self-energy on one- and two-particle properties requires to work at two-loop order, 
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which is currently inachievable. To obtain qualitative results deep inside the RC regime, 
the application of Ward identities offers an alternative. This strategy was chosen previously 
by Edwards and Hertz [14] to calculate the spin-resolved self-energy in the ordered FM 
phase. As we will show in this paper, an analogous approach can be similarly applied for 
magnetically disordered systems in the renormalized-classical regime. 

In the present paper we use the TPSC, fRG, and Ward identity approaches to get insight 
from three different points of view into the behavior of the self-energy in 2D systems at finite 
temperatures on approaching the FM phase. First, in Sect. 2 we consider the main features 
of the self-energy obtained within the TPSC approximation and compare them with the 
results of the fRG approach, which is used to study in detail the frequency dependence of 
the self-energy at van Hove (vH) band fillings. In Sect. 3 we make use of Ward identities for 
a qualitative insight into the role of self-energy and vertex corrections at low temperatures. 

II. THE SELF-ENERGY IN NON-SELF-CONSISTENT APPROACHES 

Specifically we consider the Hubbard model for N e electrons on a square lattice, 

H = ~Y1 *v4c*r + u J2 n n n n - - At ') N e (i) 

ijcr i 

where the hopping amplitude is Uj = t for nearest-neighbor sites i and j and Uj = —t' for 
next-nearest-neighbor sites (t,t' > 0); for further convenience we have shifted the chemical 
potential /i by 4t'; the corresponding electronic dispersion is 

e k = — 2t(cos k x + cos k y ) + 4t'(cos k x cos k y + 1) — /i. (2) 

Furthermore we compare some of our results for the Hubbard model with those for the 
ferromagnetic s-d model [38,39] 

H = ~ l n c \o c i<y - 1 Yl Si ' a ™ ,c \a c i°' Ji i Si " S i ~ ~ 4t ') N e ( 3 ) 

ijcr i ij 

in the weak-coupling regime < / 8t, where S« are localized-spin operators and cr^ 
are Pauli matrices; Jij = J > is the direct nearest-neighbor ferromagnetic spin exchange 
coupling. 
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A. TPSC approximation 



The TPSC approximation [9] for the Hubbard model (1) considers the dynamical spin 
susceptibility 

x(q,^) = - — - — -, r 4) 

i - ^ S pXo(q,^) 

which has the same structure as the spin susceptibility in the random-phase approximation 
(RPA), but contains an effective interaction U sp instead of the bare U. The bare susceptibility 
Xo(q,^n) is given by 

Xo(q, iu n ) = ^ J2 — ^VT7 (5) 
iv iu n — e k + e k+q 

where /k = /(£k) is the Fermi function and tu n = 2imT. The effective interaction U sp is 
determined by the sum rule [9] 

^ X(q, ^n) = n/2 - n 2 f/ sp /(4f/) (6) 

with n = N e /N being the band filling and is set to unity. The renormalization of 
the interaction U sp /U < 1 avoids the artificial divergence of the susceptibility at finite 
temperature where the denominator of Eq. (4) calculated with the bare interaction U 
vanishes. The temperature and t' dependence of U sp was extensively discussed in Refs. 
[9,26]; U sp decreases with increasing t' and decreasing temperature. 

Above the ordered ground state the inverse correlation length in TPSC 

r 1 « [1 - U spXo (QM 1/2 (7) 
monotonically decreases with temperature and becomes exponentially small, £ _1 = 

2 

C exp(— T*/T), below a certain crossover temperature T* = 4irAU sp S / (3x ), where Xo — 
Xo(Q,0), A = V 2 Xo(q, 0)| q =Q, C is a constant, and 

t xo(q,^n) 



Si = 3n/4 - 3n 2 /(8Ux ) - \ 



y - 



(8) 

T=0 



is the square of the ground-state (sublattice) magnetization per site, Q is the magnetic 
ordering wavevector, determined by the maximum of Xo(Q,0). In the following we suppose 
Q = 0, which corresponds to a ferromagnetic instability and is in particular the case for vH 
band fillings (fj, = 0) of the t-t' Hubbard model, Eq. (1), with 0.3* < f < 0.5t [32,37,26]. 
Note that we ignore here the possibility of triplet pairing, which may also arise in the vicinity 
of the FM instability [32]. 
The self-energy is given by 

UU T . 1 

E(k,ie n ) = " p V [3x(q,iw n ) - Xo(q,^n)]- —. (9) 

2M ^ lE n + lUJ n - £ k+q 

on the imaginary frequency axis [e n = (2n + l)itT are fermionic Matsubara frequencies] and 

E(k, e + i0+) = ^£ | dc[3Im X (q, c) - Im Xo (q, a;)] ^ +7o + 

on the real axis, where N B (uo) is the Bose distribution function. The factor 3 in the first 
term in the square brackets of Eqs. (9) and (10) arises from the summation over three 
(2 transverse and 1 longitudinal) spin channels, the second term in the square brackets is 
substracted to avoid double counting of the second-order diagram. 

The same expressions (9) and (10) with the replacement UU sp — > I 2 hold for the s-d model 
(3) in second order perturbation theory with respect to the electron-spin interaction /, see, 
e.g. Ref. [15]. The magnetic correlations in this case originate mainly from the exchange 
interaction between the localized spins, and the magnetic susceptibility has the same form 
as for the Heisenberg model with nearest neighbor exchange interaction J. The inverse mag- 
netic correlation length at low T is again exponentially small [40], = C exp(—JS 2 /T); 
the corresponding crossover temperature into the RC regime is T* ~ JS 2 . 

To calculate the self-energy at small e and k near the Fermi surface, we expand the 
bare susceptibility Xo & t small q and uj. For the Hubbard model at /i ^ the resulting spin 
susceptibility on the real frequency axis is given by (see e.g. Ref. [18]) 
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where the correlation length £ in TPSC is given by Eq. (7), A and B are constants which 
are proportional to U sp with a coefficient which depends on the bare spectrum e k (see Ref. 
[18] for explicit expressions). For the van Hove singularity case (/i = 0) we obtain 

xM= wte (12) 

where (cf. the T = if = result of Ref. [41]) 

r(g± ' W) = 27rtsin P2 2v9max[^,(Tt|g + g_|)V2 )t | g+g _|] ( 13 ) 

is the damping of the spin excitations, q± = q x sirup ± q y cosip, and cos2ip = It' jt. 

The form of the dynamic magnetic susceptibility of the s-d model in the paramagnetic 
phase, which is mostly determined by the local moment subsystem, is more complicated. 
However, we will see below that already the static part of the susceptibility dominates the 
spectral properties at low T < T*. In the static limit the susceptibility in the RC regime 
is expected to coincide with that determined for the AFM case within a 1/M expansion in 
the O(M) Heisenberg model [42] and has the same form as Eq. (11) with uj — 0. 

To analyze spectral properties, we first consider the inverse qp lifetime 

= -ImE(k F , i0+) = -^^N Y 3Im *^+ q )-ImXo(q, gkF+q ) 
r(k F ) 2 ^ sinh(e kF+q /T) V ; 

and the derivative of the real part of the self-energy at the Fermi level 



<9Re£(k F ,e) 



de 



UU **V I do, £[3Im X (q, u) - Im Xo (q, u>)\ N f ] + ■ (15) 



e=0 



2A^ 7 ^ " ' ' U k, c, 



The results for the leading terms at T <C i together with the results of second-order per- 
turbation theory (SOPT, which is obtained by the replacement x ~^ Xo in Eq. (14)) are 
collected in Table I, where we omit overall temperature-independent prefactors which are 
proportional to U 2 /t 2 . 
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Spectrum 




[dReZ(k F ,e)/de} £=0 




SOPT 


TPSC 


SOPT 


TPSC 


Linear (// ^ 0) 


(T 2 /t)ln(l/T) 


T£ + t£>((T/t) 2 / 3 ) 


const < 


{T/t)e + 0{{T/t)^) 


vHs (// = 0) < 


k = k vH 

Mk vll 


Tln(l/T) 
(T 2 /t)ln(l/T) 


T£ 2 + t£>((T/t) 2 / 3 ) 
T£ + t£>((T/t) 5 / 6 ) 


m(l/T) 
const < 


(T/t)e + o((T/ty/ 3 ) 

(T/t)e + 0((T/t) 2 / 3 ) 



Table I. Inverse quasiparticle lifetime T~ 1 (k.p) and [<9Re£(kp, e)/de] e=0 in second-order 
perturbation theory (SOPT) and the two-particle self-consistent approximation (TPSC) near 
the FM instability for different bare electronic spectra. 



Apparently, the imaginary part of the self-energy at the Fermi level is anomalously en- 
hanced by the correlation effects for £ 3> 1 and even tends to diverge deep inside the RC 
regime T <C T* where £ — > oo. Simultaneously, <9Re£/cte| £= o becomes positive and large. 
At energies > t^ 1 (t£~ 2 in the vHs case) the real part of the self-energy behaves as 
Aq/o;, where A = (UU sp y/ 2 S ~ T*. It is worthwhile to note that the abovementioned di- 
vergencies arise from the purely static contributions with zero bosonic Matsubara frequency 
and were previously discussed in detail for the AFM case in Refs. [6,9], where So is the 
ground-state sublattice magnetization. 

These low-energy features lead to a suppression of spectral weight in A(k, oo) = 
-(l/7r)ImS/[(cj-£ k + ReS) 2 + (ImS) 2 ] at \u\ < A (see also Ref. [26]) and to the formation 
of a two-peak structure of the spectral function with the peaks located near uj ~ ±A . 
While for the AFM case the fulfillment of the nesting condition £k F +Q = —£k F = (which 
for t' 7^ is satisfied only at the hot spots) is required, in the FM case the two-peak struc- 
ture occurs all around the FS. The difference between the FM and AFM instabilities at 
those points where £k F +Q = is evident only in the subleading terms ~ T a , where the 
exponent a depends on the dynamical exponent z [e.g. T -1 (kp) ~ T£ + t(T/i) l ~ l l z and 
<9ReE/cte| e=0 ~ (T/t)£ 2 + (T/t) 1 ' 2 ^ for the linear electronic dispersion near the Fermi level], 
z = 2 for the AFM and z = 3 for the FM case. Note that unlike Refs. [9,26] we did not 
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suppose that T <C T* in deriving the results of Table I, and in fact these results are valid 
for both, RC and QC regimes. 

The above discussed features of the self-energy keep their form away from the FS with the 
replacement e — > e — £k+Q- This holds for all |£k+Q| <C t for a linear electronic dispersion and 
for t£~ 2 <C I^ic+qI <^ t for vH band fillings. At vH band fillings and |£ic+q| ^ t£~ 2 additional 
divergent terms in the real part of the self-energy arise at the vH points k v n =(0,7r) and 
(ir, 0) with Re£(k vH , 0) ~ tk + k-£ 4 [k± = (k x —n) sm(p±k y cosy? for k vH =(vr, 0) and similarly 
for the other vH point], which flatten the bare electronic dispersion at the momenta near 
k v n, similar to the zero temperature case [45]. This flattening however is not important in 
the RC regime since the corresponding momentum region is rather narrow and ImE at the 
Fermi level is finite at finite T (and even diverges at £ — > oo). 

The dependence on e — £k+Q for t£~ 2 <C kk+cjl <C t is the origin for an important 
difference between the spectral functions near the FM and AFM instabilities away from 
the FS. In the FM case the spectral functions depend on e — e k only, implying that at 
£k — ±Ao (the condition which determines the FSs of spin-up and spin down electrons in 
the FM phase) one of the peaks of the above discussed two-peak structure is located at 
the Fermi level, making the electronic excitations at the points of the Brillouin zone with 
£k — ±Ao coherent. This indicates the existence of two "preformed" Fermi surfaces already 
in the paramagnetic (PM) phase at low temperatures T < T*. The corresponding electronic 
excitations, however, do not have any prefered spin direction and the spin symmetry remains 
necessarily unbroken at T > 0. As we show in Section III, self-consistent approaches show 
the same tendency of the spectral weight suppression at the PM Fermi surface, and the 
redistribution of spectral weight towards the energies uj ~ ±A supports therefore the 
picture described above in the non-self-consistent TPSC analysis. 
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B. fRG analysis at van Hove band fillings 



The two-particle self-consistent approximation may be insufficient close to van Hove band 
fillings since it considers only the contribution of particle-hole excitations and the other elec- 
tronic scattering channels are accounted for only by the "average" renormalized interaction 
U sp . Moreover, at van Hove band fillings the effective interaction U sp artificially tends to zero 
with decreasing temperature (see e.g. Ref. [9]). To analyze in more detail the frequency 
dependence of the self-energy in the vicinity of van Hove band fillings, we apply the fRG ap- 
proach for one particle-irreducible (1PI) functions [31] with a temperature cutoff [32]. This 
approach considers the evolution of the generating functional with decreasing temperature 
in the weak-coupling regime. The flow of the self-energy Sr(k, iu) = T _1 / 2 S(k, iuS) in the 
1PI fRG scheme is given by 

^^to S T , (16) 

where o is a short notation for the summation over momentum, frequency and spin variables 
according to standard diagrammatic rules, see e.g. Ref. [33]. The renormalization of the 
electron-electron interaction vertex Vt at one-loop order is given by 

dV T 

—— — Vt o (Gt ° S T + S T ° Gt) ° Vt- (17) 
The propagators Gt and St are defined by 

° T ^ lVn) = u^r k ] ST ^ lVn) = - 2T^^ n -e,r (18) 

The factors T 1 / 2 arise due to the rescaling of the fermion fields on removing the temperature 
dependence from the interaction term of the action [32]. Eqs. (16) and (17) have to be solved 
with the initial conditions Vr = U and S To = where T 3> t. In the non-self-consistent 
treatment in this Section we have neglected the self-energy in the denominators of the Green 
functions (18). The self-consistent RG analysis is rather involved and requires the inclusion 
of two-loop corrections which remain a challenging task for fRG techniques. 
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Since the frequency dependence of the vertices is neglected in the calculations, it is 
convenient to reinsert, following Ref. [33], the vertex from Eq. (17) into Eq. (16) to obtain 

To 



T 

To avoid the integration over temperature in the right-hand side of Eq. (19) we integrate 
by parts to obtain (cf. Ref. [33] for the momentum-cutoff scheme) 



where Tf is the final temperature where the self-energy is evaluated. Although Eq. (20) 
contains the Matsubara sums of the Green functions at different temperatures, these sums 
can be calculated by the same procedure as for equal-temperature Green functions. 

To solve Eqs. (17) and (20) numerically we divide the momentum space into 32 patches 
with the same patching scheme as in Refs. [31,32]. To calculate the self-energy on the 
real frequency axis we use analytical continuation by Pade approximants [43]. Similar to 
Ref. [11] we use the advantage of Eq. (19) that for frequency-independent vertices, after 
analytical summation over internal frequencies, the self-energy can be calculated at arbitrary 
frequencies on the imaginary axis, and therefore we choose a mesh on the frequency axis, 
which becomes denser close to iuj = [11]. 

We consider the results at the vH band filling (fx = 0) for t' = 0.45t and U — At 
where ferromagnetism is expected in the ground state [32,37]. We choose this relatively 
large value of U because the crossover temperatures T* for the FM instability are smaller 
then for the AFM case (T* < 0. It for the parameters we use), and for lower values of U 
and correspondingly lower temperatures the analytical continuation becomes increasingly 
difficult, since the size of the anomalous frequency region in the vicinity of the Fermi level 
decreases with T*. The self-energy in second-order perturbation theory (SOPT, which is 
obtained by the replacement V — > U in Eq. (19)) calculated at temperatures T = O.lt, 0.3t, 
and 0.5t for kp = (2.83, 0.07), which is the center of the first patch, closest to the (tt, 0) point, 
is shown in Fig. 2. The SOPT self-energy at other points on the FS looks similar. As in the 




(19) 



dS T 



(Gt — Gtj) ° Wt ° (Gt o St + St ° Gt) ° Vr] 



(20) 



dT 
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AFM case [11] at high temperatures T > 0.3t the self-energy has a sharp dip at the Fermi 
level due to vH singularity effects. ImE( S0PT )(kp, 0) decreases with decreasing temperature 
and ImE(kp, 0) = at T = 0. This however does not imply the validity of the qp picture 
everywhere on the Fermi surface since ImY,(k F ,e) oc e\n(l/e) for the vH points kp = (n,0) 
and (0,?r) at T = [41], although the "normal" 2D behavior ImS(k F ,e) oc e 2 \n(l/e) [44] is 
restored for other points of the Fermi surface. 

The self-energy obtained within the fRG technique for the temperatures T = O.lt and 
T = 0.3t is shown in Figs. 3-5. The results for T = 0.5t (not shown) are similar to those 
at T = 0.3t and both are close to the results of the SOPT, Fig. 2. But in contrast to 
SOPT, with decreasing temperature to T = O.lt a new sharp feature appears in a narrow 
region near the Fermi level (Fig. 4). Note that the maximum effective interaction V max = 
max{V(ki, k 2 ; k 3 , k 4 )} = 20t for this temperature. Similarly to the TPSC approach of Sect. 
Ha, the imaginary part of the self-energy has a minimum at the Fermi level instead of a 
maximum as expected for a Fermi liquid. Simultaneously, Re£(kp,c<;) has a positive slope 
near uo = 0. These pronounced anomalies in the fRG self-energy increase in size with 
decreasing T and lead to a suppression of spectral weight at the Fermi energy. The spectral 
function (Fig. 4d) has an asymmetric two-peak structure. A qualitatively similar picture is 
observed in the fourth FS patch, closest to the Brillouin zone diagonal with kp = (1.06, 0.75) 
(Fig. 5) where at low temperatures the fRG result also leads to a two-peak structure of the 
spectral function, where however the u > peak is larger than u < one. 

The magnitude of Re£(kp,0) ~ — 0.4t in the first, closest to (ir, 0) patch of the FS is 
almost the same in both, SOPT and fRG approaches, and is much larger than the cor- 
responding value found in the vicinity of the AFM instability [11]. At the same time, 
in the fourth patch of the FS, the fRG approach gives slightly positive ReE(kp,0), while 
Re£(kp,0) in SOPT approach remains almost patch-independent. Therefore, the fRG ap- 
proach, in contrast to SOPT, leads to a significant deformation of the FS near the FM 
instability at vH band fillings, which reflects a band narrowing tendency; in essence, the 
effective (renormalized) value of the next-nearest neighbor hopping t' increases towards its 
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value in the flat-band case, t'/t = 1/2 (see e.g. the discussion in Ref. [37]). This is similar 
to the earlier discussed flattening of the dispersion close to (n, 0) and (0, n) points at vH 
band fillings [45], except that in the present case the flattening affects larger parts of the 
Brillouin zone along the zone axes. The sizable value of |Re£| at the Fermi surface leads to 
an asymmetric shape of the spectral functions in the vicinity of the Fermi energy, in contrast 
to the analytical approach of Sect Ha, where this FS shift was not taken into account. 

III. THE SELF-CONSISTENT APPROACH ATT«T* 

To get qualitative insight into the role of self-consistency effects below T*, we consider 
the general form of the spin susceptibilities 

X + ~(q,iu n ) = ^ r+(k,k + q;ii/„,ii/ n + i^„)G T (k,w n )G i (k + q,iu; n + ii/ n ) 
X zz (q,iu n ) = -— ^2 F z n (k,k + q;iiy n ,iiy n + iuj n )G^(k,iiy n )G^(k + q,iu n + iiy n ) (21) 



N 



where 



r^(k,k'> n ,<) = G-\k,iv n )G-}{k\iv' n ) 

l/T l/T 

x I fe ! (- + ^')(T[^ k ,(0)cL(r)c ky (r')]>rfrrfr' (22) 







are the 3-point spin vertex functions and S°(t) = J] k cj ca (r)(T^ /3 c k+qi/ g(r). We retain the 
spin indices in this section to perform calculations in the presence of a small static external 
Zeeman field h which is set to zero at the end. To arrive at the RPA-like form of the 
susceptibilities, we follow Ref. [14] and introduce the irreducible vertices 7" CT , 

<y a (k k' IV iv' \ - r °g'( k ' k '' W "' W n) / 2 o^ 



to obtain 



* (q '^ )= l-^(q,^) (24) 
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where a, b = +, — , z, and the function 4> ab has the same structure as x ab i n Eq. (21) with 
the replacement r — > 7. Note that 

7u = 7|1 = 7ft = -7u; X + ~ = 2X 2Z 

(and similar for r) at h — 0. 

The expression for the self-energy, including the Hartree term, reads 
T1J 

s CT (k,^ n ) = -— Y^GriKK) 

k',iv' n 

TU 2 

b£<r( k » k + q; i"n, Wn + iu n )G a (k + q, iu n + iv n )x zz (q, iu n ) 
+7+_ CT (k, k + q; iv n , %v n + iu n )G_ a (k + q, + iv n )x + ~{q, iu n )\ (25) 

We neglected here the contribution of charge excitations which are not singular near the 
magnetic phase transition. 

We note that a similar expression for the self-energy can be obtained for the ferromagnetic 
s-d model (3) but in this case the susceptibilities x zz (q, iou n ) and x +_ (q, iu> n ) are determined 
by the localized-moment subsystem. This difference however will not be important in the 
following sequence of arguments and the results of this section are applicable for this model 
as well. 

As discussed in Sect. Ila, the main contribution to Eq. (22) comes from the bosonic 
Matsubara frequency iu n = 0. The vertex r" CT ,(k, k + q; iv n , iu n ) in the limit q — > [which 
we denote in the following as /(k, k; ii/ n , iv n ) and similarly for 7] is found from a Ward 
identity. The standard Ward identity (see e.g. Ref. [47]) 

lim w _ >0 r" (T /(k,k;w n ,w n + ia;) = a a aa ,[l - 9E(k, %v n + iu) / d{iu)] w=0 (26) 

is not appropriate for that purpose since it considers the opposite order of limits for the 
vertex lim^o lim g ^ T as required for the calculations in the RC regime. Instead, we use 
the identity of Ref. [14] to obtain for h — > 

r+(k, k; iv n , iu n ) = P(k, k; iv n , iu n ) = 1 + - 



h 
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h=0 



TU 2 r 



_ TC/ ^ , . , . /' 

AT 



xT(k + q, k + q; ii/ n + iu n , %v n + ^„)G 2 (k + q, icu n + w n )[x* z (q, iw n ) - X + (q, i^n)] 



+ 



7 z (k,k + q;?z/ n ,iz/ n + icj„) G(k + q, iw n + iv n )x zz (q, iu n ) \ (27) 

J ft=o J 



where we have denoted T = T + , 7 = 7+ and we have omitted spin indices. Note that in 
the third and fourth line of Eq. (27 ) we have neglected /i-derivatives of 7 + , x + ~ i an d X zz ■ 
These derivatives lead to non-singular contributions, at least for the linear dispersion case, 
which we mostly consider in the following. Using the relation (23) between T and 7 and the 
relation (21) between T and Xi we obtain from Eq. (27) the following integral equation for 
the vertex function 7 

7(k,k;u/ n ,u/ n ) ~ 1 + 
TU 2 x r 

-Jj- \ 7 ( k ' k + q ' iUn > il>n + iU ^( k + q ' k + ^ + iu) ^ ilJ ri + iuJ n) (28) 

x G 2 (k + q, iu;„ + iv n ) [x zz (q, iw n ) - X + ~ (q, iw„)] 



x (q, «w„j 



-^-7 z (k, k + q; ii/ n , ii/ n + iw n ) G(k + q, iu n + iv n ))\ 
dh J h=o j 



1 + Ux(cl, iu n ) 

In principle, Eq. (28) is the first equation in an infinite hierarchy, since the calculation of 
d^ z /dh involves d 2 ^ z /dh 2 etc. To close the system of equations in an approximate way, we 
proceed similarly to the 1/M expansion within the 0(M) generalization of the Heisenberg 
model for localized-moment systems [42,48]. Namely, we suppose that there are M — 1 
transverse spin modes and one longitudinal mode. Then the correction due to dj z /dh in 
Eq. (28) is already of order 1/M, and therefore only terms to leading order in 1/M, arising 
from the transverse spin fluctuations can be retained when calculating the derivative. To 
calculate c?7 2 (k, k; iv n , iv n )/dh we differentiate the equation for the self-energy (25) twice 
with respect to the field. Note that the terms containing d 2/ y ± /dh 2 should be retained in 
this calculation, since they are not small in 1/M and do not vanish in the q — > limit 
contrary to d^ /dh. After lengthy algebraic manipulations we finally obtain 

-^7 z (k,k;*z/ n ,^ n )l =2[l + [/ X (0,0)] 
dh J h=o 
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X nY1 7 3 (k,k + q;iz/ ri ,iz/ n + ^ n )G 3 (k + q,^ n + iz/ n )x + "(q,^ n ) +0(1/M) (29) 

Generally, the susceptibilities (24) which enter Eqs. (25) and (27) can be represented in 
a form similar to Eq. (11) or (12) at u = 0, 

R 

X + ~(q, 0) = 2x 2Z (q, 0) ~ 2 — + regular terms (30) 

where R is a constant. As in localized-moment systems with strong short-range order [46], 
the following identity 

< S ° • S '> " YM ZZ ^ °) + ^^X + -(q, 0)]e^ r = ^ exp(-r/0, (31) 
q 

which incorporates the correct long distance asymptotics of the real-space correlation func- 
tions at T <C T*, is useful to relate R and S . As a result, we obtain 

-2 M 

S 2 = -TR In £. (32) 

Since the dominant contribution to momentum sums in Eqs. (25) and (27) arises from 
the long wavelength region q < we replace 

7^'( k , k + q; iu n , iv n ) -»• 7a,<7'(k, k; u/ n , u/ n ) , (33) 

d d 
— 7 2 (k, k + q; ii/ n , ii/ n + iw n ) -> ^7 z (k, k; ii/ n , ii/ n ) . 

Equations (25), (27), and (29) together with the replacements (33) form a closed system 
for 7 and S. The remaining integrals over q are easily calculated analytically analogous to 
the non-self-consistent calculation in Section Ha. After the continuation to the real axis we 
obtain the following algebraic equations for the self-energy and the irreducible vertex 7 for 
tf _1 < \ou-e k \ <t, 

A2 7 (k,k;w,^) 



E(k,w) 



uj — £k — S(k, uj) ' 
^-2 Ag7 2 ( 

M [cj-£ k -S(k,cj)] 2 ' M[cj-e k -S(k,cj)] z 



tvv \ 1 M ~ 2 A o7 2 (k,k^,c) 2 A^ 3 (k,k^,c) 
7(k,k;w,w) = l — ~ ; . ~9. + 77~ 1 ^7T. . X14 > ( 34 ) 
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where we have introduced the ground state spin splitting A = US - Similar to the TPSC 
approach of Sect. Ha, the resulting self-energy and the vertex depend on k only through 
uj = uj — e k and are given by 

M(A 2 + uj 2 - a/cJ 2 - ai A 2 y/ZJ 2 - a 2 A 2 ) 



E(k,w) 



2 (2 + M) uj 



7(k ' k; U) = ofO ■ 1a2-2 ^ + ^ + M ) " MA 

2 (2 + M) AqCj z 

+(MAg - 2cJ 2 ) sju 2 - a.A 2 sju 2 - a 2 A 2 } , (35) 

where a lj2 = 1 +4(1 ± a/1 + M/2) /M; the branch Im-y/i > of the square roots is chosen to 
guarantee the correct analytical properties for £ and 7. On the other hand, in the absence of 
the vertex renormalization, i.e. for 7(k, k; iv n , iis n ) = 1, which is the analogue of the FLEX 
approximation [49] in our approach, we obtain for the self-energy the M-independent result 

E(k,w) = UuJ- v / cJ-2A oV / ^ + 2Ao). (36) 

To zeroth order in 1/M (i.e. at M — 00) we obtain from Eq. (35) 

A 2 A 2 
S(k, u) = ^, 7(k, k; uj, uj) = 1 - =f (37) 

i.e. the self-energy and the vertex corrections exactly cancel each other at M = 00 and the 
self-energy and the vertex are given by their second-order, non-self-consistent results. To 
demonstrate how this result changes at finite values of M, we plot in Fig. 6 the results (35) 
and (36) for M = 3 together with the corresponding spectral functions. Similarly to the 
non-self-consistent solution Eq. (37) the real part of the self-energy in the self-consistent 
solution (35) is divergent at uj — 0. This is in contrast to the solution (36), which is self- 
consistent with respect to the self-energy only, and leads to a finite positive slope of the 
real part of the self-energy at the Fermi level, dR.eT./duj = 1/2. The imaginary part of the 
self-energy has its largest absolute value at the Fermi energy and monotonically decreases 
with \uj\ for the solution (36), while it has more complicated energy dependence with a 5- 
function singularity at uj = according to the self-consistent result (35). As a consequence, 
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the solution (36), which is self-consistent with respect to the self-energy only, displays a 
one-peak structure of the spectral function, while Eq. (35) leads to a two-peak structure of 
the spectral function (Fig. 6c), which is similar to the M = oo result (37). We note that the 
gap in the spectral function at the Fermi energy in the self-consistent solution is most likely 
an artifact of the first order in 1/M result. At the same time the similarities between the 
results in zeroth- and first order in 1/M allow to conclude that the higher order corrections 
in 1/M are expected to be small and do not lead to qualitative changes. 

The one-peak structure of the spectral functions in the solution without vertex renor- 
malization is similar to that in the FLEX approximation in the AFM case, which predicts 
only a slight variation of the quasiparticle weight around the Fermi surface [7]. As discussed 
in Ref. [9], this result is a drawback of the absence of vertex corrections, and the neces- 
sity to account for frequency-dependent vertex renormalization together with the self-energy 
corrections is crucial in this case. As we have shown in this section, a similar conclusion is 
reached in the FM case; the corresponding non-trivial frequency dependence of the vertex 
is shown in Fig. 6d. 

Therefore, the solution (35) which accounts for both, self-energy and vertex corrections 
leads to the results which are qualitatively similar to those obtained in non-self-consistent 
approaches of Section II. The same arguments can be applied to the s-d model, which is 
therefore expected to lead to the similar results at T <C T*. 

IV. SUMMARY AND CONCLUSION 

We have investigated the self-energy in two dimensions in the vicinity of a FM instabil- 
ity within the TPSC, the one-loop fRG analysis, and by applying Ward identities. In all 
approaches the self-energy has a non- Fermi liquid form in narrow window \u\ < A around 
the Fermi level, and the spectral functions have a two-peak structure at low temperatures. 
The spectral weight at the Fermi energy is strongly suppressed at T < T* in both, non-self- 
consistent and self-consistent solutions, because the self-energy and the vertex corrections 
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partially cancel each other. The form of the spectral functions we have obtained in the 
paramagnetic state with strong short-range magnetic order is qualitatively similar to the 
mean-field result for the spectral functions in the FM ordered state, if the electron spin 
quantization axis is chosen perpendicular to the direction of the magnetization. Indeed, per- 
forming the simplest Stoner-like decoupling of the interaction in the Hubbard model with 
the order parameter A = (C//2)(c^c^ + c^Cjf), which corresponds to a spin alignment along 
the x-axis, we readily obtain 

« <L I c k „ >„= \ ( - - ^ 1 - - - + - - ^ J A - ;Q+ ) (38) 

where <C | 3> w denotes the Fourier transform of the retarded Green function. With the 
inclusion of a small damping 27 in the denominators the Green function Eq. (38) leads 
indeed to a spectral function, which is qualitatively similar to e.g. Fig. 6c. At the same 
time, the mean-field Green functions for the spin projection along (parallel or antiparallel) 
the magnetization axis are given by 

\ « 4 T ± 4J <*t ± c n ».= u _ £ ^ A + iQ+ • (39) 

The spectral functions in the PM phase are necessarily spin-rotation invariant. Therefore, 
we expect strong changes of the spectral functions at the zero temperature FM phase tran- 
sition only for a choice of the spin quantization axis along the direction of the ground-state 
magnetization M, while spectral functions for electrons with spin quantization axis per- 
pendicular to M change continuously with decreasing the temperature to zero (or to the 
Curie temperature T c <C t for quasi-2D systems). Note also that beyond the mean- field 
approximation, the Green functions for a spin projection antiparallel to M acquire strong 
incoherent contributions, at least for fully polarized ferromagnets [14,15]. The overall qual- 
itative picture for the evolution of the spectral properties summarizes parts of the results of 
the present paper in Fig. 7. 

The two-peak structure of the spectral functions in the vicinity of a FM instability leads 
to the formation of new coherent quasiparticles at the points of the Brillouin zone near the 
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spin-up and spin-down Fermi surfaces of the FM ordered ground state. The spin-rotation 
symmetry remains however necessarily unbroken at T > and new "preformed" qps do not 
have any preferable spin direction. This "uncertainty" of the spin direction in 2D magnets 
with strong short range order was earlier discussed in connection with NMR experiments 
[50], which are expected to give two-peak NMR spectra, which is usually characteristic for 
2-sublattice AFM order, even in the vicinity of the FM state. For quasi-two-dimensional 
systems the quasi-splitting of the Fermi surface is expected at T c < T <C A with A being 
the ground-state spin splitting. 

Translated to real space, "domains" of size ~ £ may form containing mostly electrons with 
a certain spin polarization. Unlike in localized-moments systems, however, these "domains" 
have a small electronic density and are therefore expected to be not a static but rather a 
dynamic phenomenon. Such a dynamic formation of droplets with a preferred spin direction 
near the FM instability should be distinguished from the possible phase separation into 
hole-rich and hole-poor regions [51-54] in the vicinity of an AFM instability for the almost 
half-filled band. The latter may result as a compromise to "adjust" an AFM spin structure 
to an electron density n <1 and involve charge fluctuations which are coupled to the spin 
channel for a nearly half-filled band. At the same time, the dynamic domain formation near 
a FM instability is expected solely due to FM spin fluctuations which favor the aggregation 
of electrons with certain spin orientation in regions of size ~ £. 

While the form of the spectral functions in the vicinity of the FM phase transition is 
dramatically different from the conventional quasiparticle one-peak structure, we do not 
expect sizable effects in the density of states. Indeed, since the spectral functions depend at 
low-energy on e — only, we find for \e\ <C t 



where po(e) is the bare density of states and £ m i n ,max are the energies of the bottom and the 
top of the band with respect to the Fermi level. Taking into account the two-peak structure 
of A(e) we obtain 




(40) 
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p(e) ~ ±[p Q (e - A) + p (e + A)} (41) 

i.e. the density of states is also "pre-split" at low T, but is not strongly suppressed at the 
Fermi level. 

Our results may provide a novel possibility to interpret ARPES data of layered ferro- 
magnetic materials. E.g., the described effects may be important in the interpretation of 
the results of ARPES studies of the layered manganite compound Lai+^S^-zM^Oy [3]. 
Pseudogap structures are observed in this material both, above and below Tq- It was sug- 
gested that these structures are possibly produced by an accompanying charge order [55], 
or by phase separation [56]. The FM fluctuations however might be responsible for the part 
of the pronounced shift (~250 meV) of the spectral weight maxima away from the Fermi 
energy at the points of the Brillouin zone near the FM Fermi surface with T above the Curie 
temperature T C - 

The described effects may be also important for the normal state of some unconventional 
superconductors, where ferromagnetic fluctuations are expected to be important (c.f. UGe 2 , 
Sr 2 Ru04). UGe2 has long range magnetic order in the ground state and therefore may 
also show a quasi-splitting of the FS above T c . Although long range FM order is absent in 
strontium ruthenate Sr 2 Ru0 4 , FM order is induced by a small amount of La doping [57]. 

Taking into account the strong self-energy and vertex renormalization in the vicinity of 
the FM phase transition, the criteria for ferromagnetism [18] should be reconsidered. This 
especially concerns quasi-2D systems, where these renormalizations are expected to be the 
strongest. However, even in 3D systems, where the finite-temperature divergencies of the 
self-energy at the Fermi level are only logarithmic in £ (instead of a power law divergence in 
2D), they can also lead to additional renormalizations of the Stoner criterion in comparison 
to those considered within the paramagnon (spin fluctuation) theory of Refs. [18,19]. 

So far unresolved is the behavior of the self-energy in the quantum-critical regime, if it 
exists at all for the FM phase transition due to the presence of non-analytic contributions. 
This behavior depends crucially on the temperature dependence of the correlation length; the 
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£ 2/3 dependence of the self-energy at T = however already implies a non-trivial temperature 
dependence of the quasiparticle scattering rate at finite temperatures. The simple ansatz 
for the magnetic susceptibility we have used in the Ward identity approach is not justified 
in this case, and an alternative approach has yet to be developed. 
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FIGURES 




FIG. 1. Schematic picture for the three different temperature regimes, as proposed in Rcf. 
[25]: (1) the renormalized classical (RC) regime above the ordered ground state (indicated by 
the bold line, T* is the corresponding crossover temperature scale discussed in the text), (2) the 
quantum disordered (QD) regime above a disordered ground state, and (3) the quantum critical 
(QC) regime. The arrow shows the direction, in which the evolution of the spectral properties is 
traced in this paper; "fRG" and "Ward identities" in the figure mark the intermediate (T > T*) 
and low temperature (T <C T*) regimes, where the corresponding approaches are applied. The 
two-particle self-consistent approach is applied in both regimes. 
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FIG. 2. The real (a) and imaginary (b) parts of the self-energy in second-order perturbation 
theory (SOPT) at t'/t = 0.45, U = 4i, fj, = (vH band filling), and different temperatures. 
Arrows mark the energy of the lower- (e m i n = —4 + 8t') and upper- (e max = 4 + 8t') edges of the 
noninteracting band, the cross marks the energy, opposite to the upper edge of the band, — e max . 
The inset shows the behavior at low energies for T = O.lt. 
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FIG. 3. The functional renormalization group results for the self-energy (a-c) and the spectral 
function (d) in the first Fermi surface patch, closest to the (tt, 0) point at t' /t = 0.45, U = 4t, 
fi = 0, and T = 0.3t. 
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FIG. 4. Same as Fig. 3 for T = O.lt. 
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FIG. 5. Same as Fig. 3 for T = O.lt in the fourth Fermi surface patch closest to the Brillouin 
zone diagonal. 
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FIG. 6. The real and imaginary parts of the self-energy (a,b), the spectral function (c), and the 
vertex function 7 (d) at T <C T* in the fully self-consistent Ward identity approach (solid lines) 
and the FLEX-like approach, which is self-consistent with respect to the self-energy only (dashed 
lines) ftS ft function of lj = u — £k- 
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FIG. 7. Schematic picture for the evolution of the spectral function near the ferromagnetic 
phase transition. The spectral functions at T = are shown at the paramagnetic Fermi surface. 
The spectral functions at the Fermi surfaces of spin-up and spin-down electrons are expected to 
shift by ±A with respect to those of the paramagnetic Fermi surface, as indicated by arrows. 
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